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We present a finite size scaling analysis of Monte Carlo simulation results on a four dimensional 
' Ising spin glass. We study chaos with both coupling and temperature perturbations, and find the 

, same chaos exponent in each case. Chaos is investigated both at the critical temperature and below 

^\ • where it seems to be more efficient (larger exponent). Dimension four seems to be above the critical 

7—i ' dimension where chaos with temperature is no more present in the critical region. Our results are 

,—1 , consistent with the Gaussian and bimodal coupling distributions being in the same universality class. 

^ ■ 

^ '■ I. INTRODUCTION 

T— I ■ 

I— I ■ While the nature of the spin glass phase is still controversial, there is a property, namely static chaos, that emerges as 
Ph ' a comraon feature ai different spin glass models. This chaotic behavior has been studied within mean- field theorytl u, 
5^ ' scalingQ or dropletB theories and using a real space renormalisation-group approachQ'Ll. Chaos means that the frozen 
do random equilibrium state of the spin glass phase is completely reorganised by a small change in an external parameter, 
such as temperature (T) or magnetic field. Chaos with respect to a slight change in the the couplings, a so-called 
random perturbation, has also been studied. In the case of a temperature preturbation, for example, the spin-spin 
, correlation varies chaotically with T, the larger the distance between the spins, the bigger this effect. The temperature 
scale at which the correlation function varies is 



. ATl- L^^, (1) 

O ' which defines the Lyapunov or chaos exponent ^. The existence of such a scaling and of C has been shown within the 
models quoted above and a summary of the findings will be presented later in this introduction. Here we study chaos 
by Monte Carlo (MC) simulations of a four-dimensional spin glass, and we consider the Edwards- Anderson model 
~^ ; with Gaussian coupling distribution. Some of the questions adressed in this paper are: 

^ • Arc all perturbations equivalent ? Chaos with two different random perturbations and a temperature change 

, have been studied at the critical temperature and are found to give a unique exponent, within the uncertainties. The 
amplitudes of the effect are not the same, however, and a temperature perturbation is more difficult to see numerically 
than the random perturbations. 
■ • Are there different chaos exponents at the critical temperature and in the spin glass phase ? We get rather 
f — ■ ^ different values of C, at Tc and below. This seems to support previous findings that chaos is more effective in the spin 
0^ . glass phase (larger exponent )I3. 

' A brief reminder is needed to compare these results with ones o£ other models. Chaos in spin glass phase was 
^ , first pointed out in mean-field theory by, ParisitJ, and later confirmed]. A loop-expansion around Parisi's solution for 
^ ■ dimensions d > 8 allowed Kondor et al.0 to show chaos with magnetic field and also with temperature. Two distinct 
_L , exponents ^ resulted while the temperature perturbation has a smaller effect. l-. 
^ ■ For low dimensional systems, chaos has been extensively studied within the scaling theory of Bray and MooreQ and 
Q ' the droplet theory of Fisher and HuseQ which are based on a real space renormalisation group approach. The latter 
O ' permits one to determine chaos exponents for arbitrary perturbations and in various points of the phase diagram. 
^ ' Chaos has been shown to be characteristic of each fixed point of the diagramQ. In particular, one gets different 
• ^ . exponents at the critical temperature (Tc) (between spin glass and paramagnetic phases) and in the ordered phase 
' (T < Tc), in dimension thrge. One can argue that there is a critical dimension above which there is no chaos with 
^ I temperature anymore at TcU. In the Migdal-Kadanoff (MK) framework which is the renormalisation scheme that has 
■ - - ' been mainly used, dimension four is above this critical dimension. This leads to an additional question for our study: 
• What is this critical dimension above which chaos with temperature no longer exists at the critical temperature? 
Our Monte Carlo simulation indicates that four is close to but higher than this critical dimension. 

Moreover, within MK, the exponents from a temperature perturbation or a random perturbation have to be the 
same, keeping everything else fixed, (in agreement with this paper's results). Finally, a magnetic field perturbation 
has been less studied but seems to lead to a different exponenta. As it was the case with mean-field approach, chaos 
with temperature is less effective (smaller exponent) than with magnetic field. A major difference, however, is that a 
magnetic field destroys the spin glass phase and chaos is studied in the paramagnetic phase close to zero field while, 
in mean field, one stays below a non-zero critical field. 
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MC simulations of the Edwards- Anderson model allow comparisons with both previous approachesB. Ritort's 
simulationalj of the mean-field version of this model confirms the analytical results of Kondor. In dimension two, 
where Tc is zerOj-we have been able to observe chaos with both temperature and random perturbation and found a 
unique exponentHJ. We are now interested in higher dinaeiisions where there is a spin glass phase. Since dimension 
four has a clear transition in contrast to dimension thrcdlil, and also because there has already been some studies of 
chaos for ±J random couplingalj, we focus on a Gaussian coupling four-dimensional system. Thus we are able to 
answer to the question : 

• Are exponents the same for different coupling distributions ? The results are consistent with the Gaussian (this 
paper) and bimodal (Ritort et al.O) distributions lying in the same universality class. _ _ 

Chaos with temperature or magnetic field are believed to explain cycling temperaturecJ or marantic fielctj ex- 
periments on spin glasses. Such an experiment has been done recently on a disordered ferromagnettJ. This can be 
understood whithin renormalisation group approachQ where chaos is shown to be present when the coupling distri- 
bution is shifted towards ferromagnetic couplings. It could be interesting to check that chaos can also be seen with 
Monte Carlo simulations. 



II. THE MODEL 



The Hamiltonian (for Ising spins {St] and nearest neighbor couplings { Jij}), is 

H = — ^ JijSiSj , (2) 

where the couplings are drawn from a Gaussian distribution with zero mean and variance Jf^ equal to unity. The s p in s 
lie on a four-dimensional cubic lattice of linear size L with periodic boundary conditions. As in previous studiestil o, 
the basic quantity is the replica overlap between two copies (replicas a and b) of the system 

1 

1ab = <S'|°''S'f ^ . (3) 
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From this the Binder ratio is computed, 



'=2 



4\ 1 
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(4) 



where (• ■ •) denotes both the average over disorder and the statistical mechanics (Monte Carlo) average. It is dimen- 
sionless and has a finite size scaling which allows for calculation of the critical temperature 

ii/'^(T-r,)V (5) 



9 = 9 

A similar approach has been introduced by Ritortlll to investigate chaos with random perturbation. One has now 
to compare two copies with correlated coupling sets. We obtain the perturbed couplings {Jij} from the unperturbed 
ones {Jij} in two ways. First, 

(1) 4 = 'IilI^M^ , (6) 

where is a Gaussian random variable with zero mean and unit variance. Secondly, we consider the perturbation 
that changes the sign of a small fraction of the couplings, 

(2) J-j = — Jij with probability p (7) 

and J-j = Jij otherwise. With ± J coupplings, one can only study case (2)0. In order to compare both perturbations, 
one can determine how the two sets of couplings are correlated. One gets 

AJ2 



(1) J^J•n, 



(2) J,,4. = l~2p. (8) 



2 



where (1) is expanded for small A J. r-il— i 
A measure of the spin reorganisation under small perturbation is the chaos parametedlj'cj 

where q.j.j' is given by Eq.(H) with two copies a and b having now slightly different couplings. Note that qjj = q.u' 
since {J'ij}, and the {Jij} have the same distribution. The scaling of this quantity, at fixed temperature, leads to a 
chaos exponent for each random perturbation (1) and (2), 

r^j = f^j (L^AJ) (10) 

For a temperature perturbation, one can follow the same pathwayQ and define a parameter r^j, like in Eq. (||) 



(11) 



where qT^T_ is given by (^ with two copies a and b having slightly different temperatures, T_ and T+, which are 
equally shifted from a reference temperature T, that is, 

T±^T±AT/2 . (12) 

Finally, an analoguous finite size ansatz is applied to this parameter 

(L^AT) , (13) 

which could lead to another chaos exponent C. (although we use the same symbol). 



III. RESULTS AND DISCUSSION 



A. Chaos at 

A plot of g as a function of T for various sizes L shows an intersection at = 1.8, see Fig.l. Kapwing this permits 
one to determine the critical exponent v from a plot of g as a function of (T — Tc)L^^'^ , as usualll^O. In such a plot, 
see Fig. 2 , a collapse of all data to a single curve is obtained for 

Tc = 1.8 ±0.05, 1/ = 0.87 ±0.15 . (14) 

This is in agreement with previous MC simulationJ§'lll. 

Scaling plots of the chaos parameter for the three perturbations are shown in Fig. 3, 4 and 5 where the reference 
temperature is precisely Tc. In the case with random perturbations, Eqs. (|6[ |^), and temperature change given by 
Eq. (|l|) with T = Tc, one gets 

chaos with AJ : C = 0.85 ±0.1 

chaos with AT : C = 0.95 ± 0.2 . (15) 

One can check that exponents from both random perturbations are equal (see Fig. 3 and 4). The data for chaos with 
T in Fig. 5 has larger statitical errors and also deviates less from unity; in other wordSp-|the amplitude of chaos with 
T is smaller which is also the case with other simulationstHi and mean-field calculationscl. 

We conclude that the exponents for chaos with A J and AT at Tc, given by Eq. (|l5|), are equal within the error 
bars. We use this to compare results from several sources for chaos with AJ and AT together in Table I. 

In order to define chaos with AT in the critical region it is necessary that C > 1/i^ so that the typical temperature 
interval on which the spin correlation varies, AT^ ~ T~^, is smaller than the critical temperature range given by 
T — Tc| ~ L~^/^. Both exponents depends on dimension. In our four dimensional system \/v = 1.15, and the 
inequality does not seem to be satisfied (or there might be an equality) ruling out chaos with T in the critical region. 
This is also the case with MKRG which gives 1/v = 0.68 (to be compared with C, = 0.53 in Table I). According to 
MKRG, chaos at the critical temperature is present in dimension three, \/v ~ 0.35; this is not yet confirmed by 
Monte Carlo simulations. 

Finally, one can see from Table I, in dimension four (numbers in brakets), that the exponents for Gaussian and ± J 
distributions are equal (error bars are of order ±0.1). 
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TABLE I. Values of the chaos exponent is shown for various dimensions. In some cases we have ( both in the spin glass 
phase, T <Tc, and in the critical region, at Tc (numbers in brakets). The models are Migdal Kadanoff Renormalisation group 
calculations (MKRG), mean-field expansions around Parisi's solution, Monte Carlo simulatie*is for Gaussian (G) coupling 
distributions (this paper's results in dimension four, and our previous results in dimension twctj), and bimodal (±J) coupfing 
distribution (Note that Ritort et al. used another symbol X = 2/() and finally exact ground state numerical calculations. The 
chaos exponent has been determined both with random and temperature perturbations except when it is mentioned A J or AT. 



dimension 
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T, = 


T < T,(r,) 


T < T,{T,) 


T < r. 


MKRtX 
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(AT) 








1.0 


Monte CarloE 


(±J) 
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1.2 (0.85) 




Ground Statet. 


3 (A.T) 


0.95 

















B. Chaos in the spin glass phase 



We also calculate the chaos exponent in the spin glass phase. From a previous numerical result of Young et al. 




temperature T = 1.4 and sizes L > 3 seems to probe the ordered low temperature phase. The data for T — 1.4 is 
shown in Fig. 6 from which we estimate 

chaos with A J, case (1) : C = 1.2 ±0.1. (16) 

This exponent is larger (although we cannot rule out an equality) than the one at T^- This indicates that chaos is 
more efficient at low temperature in agreement with other sources quoted in Table I, see the dimension four column. 



C. Conclusions 



We find that the Gaussian and bimodal distrihuitions seems to be in the same universality class (in contrast to 
numerical calculations of other critical cxponentsEJ). p- ■ 

Our results on chaos in four dimensions (this papec) and in two dimensions^ are in qualitative agreement with 
those of the real-space-renormalisation-group approachQ. Some of the results obtained in this approachLl apd observed 
numerically are the following. First, chaos is also present in the paramagnetic phase (e.g., in dimension twoEJ). Second, 
a temperature perturbation generates a random perturbation and thus gives the same chaos exponent (confirmed 
numerically in both two and four dimensions). Moreover, chaos with temperature occurs inside the critical region 
if the dimension is smaller than a limit dimensionEI, called d+, with 3 < d+ < 4 (in agreement with this paper's 
results). The chaos exponent is smaller at Tc than in the spin glass phase (see Table I). More generally, it has been 
arguedQ that a chaos exponent can be assigned to each fixed point of each random system, and can appear in physical 
quantities when a perturbatian of the thermodynamic parameters couples to a random perturbation. The disordered 
ferromagnet is one exampletZI, and it could be interesting to perform a Monte Carlo simulation of a spin glass with 
random ferromagnetic couplings. 
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FIG. 1. Data for g against T. Prom the intersection, Tc is estimated to be 1.8. 
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. 2. A scaling plot of the data in Fig.l scaled according to Eq. which gives u — 0.87 {l/f ~ 1.15). 
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FIG. 3. A scaling plot of r&j with random perturbation, case (1), at Tc- The perturbation lies in the range 0.05 — 0.4. 
Trying different values, our best estimate is C = 0.85. 
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FIG. 4. A scaling plot of taj with random perturbation, case (2), at Tc- To compare these data with Fig. 3, we use Eq. 
(^) and define A J = 2p^^^ . The perturbation p lies in the range 0.001 — 0.1 which is A J = 0.063 — 0.63. The chaos exponent 
is C = 0.85. 
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FIG. 5. A scaling plot of tat at Tc. The perturbation, AT, lies in the range 0.05 — 0.25. The chaos exponent is C = 10. 
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FIG. 6. A scaling plot of Taj with random perturbation, case (1), at T = 1.4 in the spin glass phase. The perturbation lies 
in the range 0.05 — 0.4. The chaos exponent is C = 1-2. 
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